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ABSTRACT 

Context. X-ray observations have unveiled the existence of enigmatic point-like sources at the center of young (a few kyrs) 
supernova remnants. These sources, known as Central Compact Objects (CCOs), are thought to be neutron stars produced by 
the supernova explosion, although their X-ray phenomenology makes them markedly different from all the other young neutron 
stars discovered so far. 



> 

' Aims. The aim of this work is to search for the optical/IR counterpart of the Vela Junior CCO and to understand the nature 
, of the associated Ha nebula discovered by Pellizzoni et al. (2002). 

■ Methods. We have used deep optical {R band) and IR (J, H, Ks bands) observations recently performed by our group with 
the ESO VLT to obtain the first deep, high resolution images of the field with the goal of resolving the nebula structure and 
pinpointing a point-like source possibly associated with the neutron star. 

Results. Our _R-band image shows that both the nebula's flux and its structure are very similar to the Ha ones, suggesting that 
the nebula spectrum is dominated by pure Ha line emission. However, the nebula is not detected in our IR observations, whick 
• • , makes it impossible to to constrain its spectrum. A faint point-like object (J > 22.6, H ~ 21.6, Ks ~ 21.4) compatible with the 
^ ' neutron star's Chandra, X-ray position is detected in our IR images {H and Ks) but not in the optical one [R ^ 25.6), where 
k>( [ it is buried by the nebula background. 
j_j . Conclusions. The nebula is most likely a bow-shock produced by the neutron star motion through the ISM or, alternatively, a 
' photo-ionization nebula powered by UV radiation from a hot neutron star. A synchrotron nebula, powered by the relativistic 
particle wind from the neutron star, is the less likely interpretation because of its non-detection in the X-rays and of the apparent 
lack of continuum emission. The candidate CCO counterpart could be the neutron star itself, a fallback disk left over from the 
supernova explosion, or a mid M-type star, or later, at a distance of 2 kpc. 

Key words. Stars: neutron. Stars: individual: CXO J085201.4-461753 
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1. Introduction 

Ground breaking observations with the current generation 
of X-ray observatories have recently unveiled unexpected 
and puzzling aspects of the phenomenology of Isolated 
Neutron Stars (INSs). It has become clear that not all 
INSs are born as fast spinning, active radio pulsars. A 
handful of isolated compact sources, most probably young 
INSs, do not show any "standard" pulsar activity. They 
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are known as "Central Compact Objects" (CCOs) in su- 
pernova remnants (Pavlov et al. 2004), enigmatic com- 
pact sources lying very close to the center of young (0.3- 
20 kyrs) supernova remnants, showing only thermal-like, 
often unpulsed X-ray emission. They are supposed to be 
the youngest members of the radio-quiet NS family which 
also includes the Anomalous X-ray Pulsars (AXPs) and 
the Soft Gamma-ray Repeaters (SGRs), which make up 
the class of the magnetars (Woods & Thompson 2006), 
the INSs with X-ray dim thermal emission (XDINSs, 
Haberl 2005), as well as high-energy pulsars like Geminga 
(Caraveo et al. 1996). At present, we cannot even exclude 
that such sources represent "the rule" among INSs, radio 
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pulsars being "the exception" , but much easier to be ob- 
served. Their overall properties are largely unknown. First 
of all, the lack of observable radio emission has still to be 
explained. It could be possibly due to an unfavourable 
beaming, to inhibition by an ultra-high magnetic field, as 
invoked for the magnetars, or to quenching by accretion 
from a residual disk. The nature of the CCO thermal X- 
ray emission is also unclear. It could be produced by the 
cooling NS surface, by hot polar caps, by the NS magnetic 
field decay, or by accretion from a residual disk. Their re- 
lationship (if any) with other classes of radio-quiet INSs 
is also debated. 

As in the case of other radio-quiet INSs, optical/IR 
observations could provide a major step forward in the 
understanding of the nature of the CCOs. Sanwal et al. 
(2002) identified a field M dwarf detected by EST as a 
likely IR counterpart to the RCW103 CCO, possibly mak- 
ing it the first X-ray Binary associated to a SNR. The 
identification, however, is still unconfirmed. IR observa- 
tions of other CCOs have been recently carried out (Fesen 
et al. 2006; Wang et al. 2007), but no viable candidate 
counterpart has been identified. An IR counterpart to the 
PKS 1209-52 CCO was proposed by Pavlov et al. (2004) 
but the identification has been recently discarded both by 
Mignani et al. (2007a) and Wang et al. (2007) on the basis 
of its significant offset with respect to the Chandra posi- 
tion. The unsuccessful identification scoreis explained by 
the fact that most CCOs are rather distant and heavily 
absorbed, and therefore they are very difficult targets at 
optical/IR wavelengths. 

The only possible exception is the X-ray source at 
the center of the RXJ 0852.0-4622 supernova remnant, 
indeed one of the least studied CCOs. RXJ 0852.0-4622 
(G266.1-1.2) is a very young (a few thousands years) shell- 
like supernova remnant discovered in the ROSAT All Sky 
Survey (Aschenbach 1998), whose projected position is co- 
incident with the southeast edge of the more extended 
Vela supernova remnant, hence dubbed "Vela Junior". 
However, its estimated distance of ~1 kpc (Slane et al. 
2001) puts it beyond the Vela supernova remnant and 
thus it rules out any intriguing association between the 
two. The CCO in Vela Jr. (AX J0851.9-4617.4) was dis- 
covered by ASCA (Slane et al. 2001) and studied with 
BeppoSax (Mereghetti 2001) and later with Chandra, 
which also provided its sub-arcsec position (Pavlov et al. 
2001). The X-ray emission of the Vela Jr. CCO, hereafter 
CXO J085201. 4-461753, is characterized by a thermal-like 
spectrum with kT 404 eV (A'^h ~ 3.510^^ cm"^ ) and an 
emitting radius of « 0.28 km (at 1 kpc), with no evidence 
of pulsations (Kargaltsev et al. 2002). First optical inves- 
tigations of the CXO J085201.4-461753 field were carried 
out by Pavlov et al. (2001) with the CTIO 0.9m telescope, 
and by Mereghetti et al. (2002) using archival observa- 
tions obtained with the ESO/MPG 2.2m telescope. The 
derived optical flux upper limits of B 23 and R ^ 22.5 
(Mereghetti et al. 2002) imply a very high X-ray to op- 
tical flux ratio which virtually certified CXO J085201.4- 
461753 as the compact remnant of the Vela Jr. supernova 



explosion. Although no optical candidate counterpart was 
detected by Mereghetti et al. (2002), a short expo- 
sure unveiled the presence of an extended emission feature 
(~ 6" diameter) positionally compatible with the CXO 
J085201. 4-461753 error circle. The reality of this feature 
was confirmed soon after by Pellizzoni et al. (2002) us- 
ing a deeper Ha exposure (^ 3 hours) taken as a part of 
the Southern Galactic Plane and Magellanic Clouds dig- 
itized plate survey (Parker & Phillips 1998) carried out 
with the AAO UK Schmidt Telescope. Unfortunately, the 
low spatial resolution of the Ha image (1" /pixel) did not 
allow to resolve the feature structure and to favour one 
of the two alternative proposed interpretations i.e., either 
a velocity-driven bow-shock nebula, or a photo-ionization 
nebula, powered by UV radiation from a hot INS, and 
therefore to better constrain the properties of the CCO. 
Here we present the results of the first deep optical and IR 
observations of the Vela Jr. CCO recently performed with 
the ESO VLT. Optical and IR observations are presented 
in §2 and §3, respectively while the results are discussed 
in §4. 

2. Optical Observations 

2.1. Observations Description 

We have performed deep optical observations of the Vela 
Jr. CCO with FORSl (FOcal Reducer Spectrograph), a 
multi-mode instrument for imaging and long-slit/multi- 
object spectroscopy mounted at the second Unit Telescope 
(UT2) of the VLT (Paranal Observatory). The obser- 
vations were carried out in Service Mode on January 
15th and 17th 2005 through the R Bessel filter (A = 
6570 1; AA = 1500 A). The four port read out mode and 
high gain were choosen as the default instrument config- 
uration. To minimize the light pollution from "star Z" of 
Pavlov et al. (2001), located ~ 1'.'5 away from our tar- 
get, we split the integration time in 20 exposures of 260 s 
each. In order to achieve the best possible spatial resolu- 
tion, useful to disantangle any point source contribution 
from the PSF wings of star Z, FORSl was used in its 
High Resolution (HR) mode (1x1 binning), with a corre- 
sponding pixel size of O'.'l. The very bright stars HD 76060 
{V = 7.88) and Wray 16-30 {V = 13.8), respectively w 35" 
and « 25" away from our target (Mereghetti 2001), were 
masked using the FORSl occulting bars. The observations 
were collected with an average seeing of ~ 0'.'8 — 0'.'9 and 
an airmass ~ 1.3 (see Table [T]) on both nights and with a 
fractional lunar illumination of 57% and 47 % for the first 
and second night, respectively. 

2.2. Data Reduction and Analysis 

Usual reduction steps to remove instrumental signatures 
were applied through the ESO FORSl data reduction 
pipeline^ using bias frames and twilight flat fields. Then, 

1 www.eso.org/observing/dfo/quality/FORSl/pipeline 
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Table 1. Summary of the VLTFORSl i?-band observa- 
tions of the Vela Jr. field with the number of exposures 
per night, the single exposure time, the average seeing and 
airmass values. 



yyyy.mm.dd 


N 


Exp. (s) 


Seeing 


Airm. 


2005.01.15 


10 


259 


0'.'83 


1.3 


2005.01.17 


10 


259 


0'.'93 


1.3 



single reduced images were combined to filter out cosmic 
rays. Since no dither was applied between single exposures 
to preserve the masking of the bright stars in the field, no 
image registration was required before the frame stack- 
ing. The photometric calibration was performed through 
the observation of standard stars from the field PG1323-86 
(Landolt 1992), yielding i?-band zero points of 26.86±0.05 
and 27.48 ±0.03 for the first and the second night, respec- 
tively. We note that the first night was not photometric, 
with the computed zeropoint deviating by ^ 0.6 magni- 
tudes from the tren Thus, we have re-normalized our 
photometry using as a reference the zeropoint of the sec- 
ond night which, instead, is very well consistent with the 
zeropoint trend. 

The coordinates of the Vela Jr. CCO were measured with 
a good nominal accuracy (0'.'6) by Pavlov et al. (2001) 
who report: a( J2000)=08'' 52™ 01.38", (5(J2000)== -46° 
17' 53'.'34. We have re-analyzed the same Chandra obser- 
vation and we have obtained: a( J2000)=08'' 52™ 01.37", 
(5( J2000)= -46° 17' 53'.'50, i.e. with a shght offset south- 
west but statistically consistent with the previous ones. 

The astrometry on the FORSl image was computed 
using as a reference 26 stars selected from the 2MASS 
catalogue. The pixel coordinates of these stars (all non 
saturated and evenly distributed in the field) were mea- 
sured by gaussian fitting their intensity profiles by using 
the specific function of the GAIA (Graphical Astronomy 
and Image Analysis) toolf] while the fit to the RA, DEC 
reference frame was performed using the Starlink package 
ASTROlvifl. The rms of the astrometric solution turned 
out to be « 0'.'12 per coordinate. After accounting for the 
0'.'2 average astrometric accuracy of the 2MASS, the over- 
all uncertainty to be attached to the position of our target 
is finally 0'.'65. 

A 8" X 8" zoom of the FORSl i?-band image of the 
field is shown in Fig. [T]-left. No point-like source appears 
at the Chandra position, thus setting an upper limit of 
R ~ 25.6 (3cr) on the CCO counterpart i.e., about 3 mag- 
nitudes deeper than the one obtained by Mereghetti et al. 
(2002). Instead, a compact optical nebula is cleary de- 
tected. We exclude that this nebula is an artifact due 
to a PSF anomaly of star Z, to a defect in the image 
flat fielding, to ghost images, or to any other instrumen- 
tal effect. Both its position and extent though (Fig. [T]- 

www.eso.org/ observing/ dfo / quality/FORS 1/qc /zeropoints 
^ star-www.dur.ac.uk/ pdraper/gaia/gaia.html 
* http:/ /star-www. rl.ac.uk/Software/software. htm 



Table 2. Summary of the VLTNACO J, H, X-band ob- 
servations of the Vela Jr. field with the number of expo- 
sures per filter (NxNDIT), the Detector Integration Time 
(DIT) , the average seeing and airmass values. 



yyyy.mm.dd Fil. NxNDIT DIT (s) Seeing Airm. 



2006.05.23 




5x19 


24 


0'.'78 


1.18 




H 


2x19 


60 


O'.'Sl 


1.32 


2006.05.24 


Ks 


2x20 


60 


0'.'74 


1.19 




J 


2x21 


60 


0'.'63 


1.25 



right), are consistent with the one of the Ha nebula seen 
by Pellizzoni et al. (2002), which clearly indicates that 
they are the same object. The nebula brightness in the R- 
band is ~ 22.8 magnitude arcsec"^ in the central part. 
Unfortunately, no color information is available on the 
nebula and, thus, it is not possible to constrain its spec- 
trum. However, we note that the i?-band flux of the nebula 
is (3.7±0.3) X 10"^'' erg cm~^ s~^, which is consistent with 
the value of (3.0 ± 0.6) x 10^^^ erg cm~^ s~^ measured 
in Ha by Pellizzoni et al. (2002). This suggests that the 
nebula spectrum has, at least in the R band, no significant 
continuum components and is dominated by pure Ha line 
emission. 

3. IR Observations 

3.1. Observations Description 

Deep IR observations of the Vela Jr. CCO were performed 
in visitor mode on May 23rd and 24th 2006 with NAos 
COnica {NACO), an adaptive optics (AO) imager and 
spectrometer mounted at the fourth Unit Telescope (UT4) 
of the VLT. In order to provide the best combination be- 
tween angular resolution and sensitivity, NACO was op- 
erated with the S27 camera with a corresponding field of 
view of 28" x 28" and a pixel scale of 0"027. As a ref- 
erence for the AO correction we have used the GSC-2 
star S1331311130291 {V = 15.3), located 11"3 away from 
our target. The Visual {VIS) dichroic element and wave- 
front sensor (4500 — lOOOOA) were used. Observations were 
performed in the J(A = 12650 A; AA = 2500 A), H{X = 
166001; AA = 3300 A) and K,{X = 21800 A; AA = 3500 A) 
filters, with a total net integration time of about 2300 s per 
band. To allow for subtraction of the variable IR sky back- 
ground, each integration was split in sequences of short 
randomly dithered exposures with Detector Integration 
Times (DIT) of 24 and 60 s and NDIT repetitions along 
each point of the dithering pattern (see Table [2|) . The in- 
strument readout mode was selected according to the DIT 
in order to minimize the read out noise. For all observa- 
tions, the seeing conditions were on average below ~ 0"8 
and the airmass was better than 1.3, allowing for an op- 
timal use of the NACO adaptive optics. Unfortunately, 
in the first night the AO correction degraded signifi- 
cantly towards the end of the iJ-band exposure sequence. 
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Fig. 1. (left) Inner section (8" x 8") of the VLT/FORSl i?-band image of the Vela Jr. CCO. North to the top, East to 
the left. The circle (0'.'7 radius, la) corresponds to the recomputed Chandra position uncertainty after accounting for 
the accuracy of our astrometric solution (see §2.2). Star Z of Pavlov et al. (2001) is labelled. The cross indicates the 
position of the candidate CCO counterpart (see §3.3). (right) Isophotal contours from the UKST Ha image of Pellizzoni 
et al. (2002) overplotted. Dark and light isophotes have a linear and geometric spacing in intensity, respectively. The 
coincidence of the diffuse feature both in position and in morphology is evident. 



Night (twilight flat flelds) and day time calibration frames 
(darks, lamp flat fields) were taken daily as part of the 
NACO calibration plan. Standard stars from the Persson 
et al. (1998) fields were observed in both nights for pho- 
tometric calibration. 



3.2. Data Reduction and Analysis 

The data have been processed using the ESO NACO 
pipelined and the science images reduced with the pro- 
duced master dark and flat field frames. For each band, 
single reduced science exposures have been then combined 
by the NACO pipeline to produce cosmic-ray free and 
sky-subtracted images. The photometric calibration 
pipeline yielded average zero points of 23.87 ± 0.07 (J), 
23.81 ±0.05 (iJ), 22.93 ±0.03 (if^) for the first night and 
23.98 ± 0.04 (J), 23.85 ± 0.04 (iJ), 22.85 ± 0.03 (K,) for 
the second night. 

The astrometric calibration of the NACO images was 
performed using the same procedure described in §2.2. 
However, since only four 2MASS stars are identified in 
the narrow NACO S27 field of view, we have computed 
our astrometry by using as a reference a set of ten sec- 
ondary stars found in common with our FORSl _R-band 
image, calibrated using 2MASS (see §2.2). The rms of the 
astrometric solution turned then to be « 0'.'06 per coordi- 
nate. By adding in quadrature the rms of the astrometric 
solution of the FORSl image (« 0'.'12 per coordinate) and 



www.eso.org/observing/dfo/quality/NACO/pipeline 



after accounting for the astrometric accuracy of 2MASS 
we thus end up with an overall accuracy of 0'.'37 per co- 
ordinate on the NACO images astrometry. Assuming the 
same Chandra positional error as before (0'.'6), the overall 
uncertainty on the target position on the NACO images 
is thus 0'.'7. 



3.3. Results 

We have identified a very faint, apparently point-like, ob- 
ject close to the southeast edge of the Chandra CCO error 
circle (see Fig. [2]). Although it is detected with a marginal 
significance (^ 4cr), we exclude that this object is a fake 
detection due to a background fluctuation, to a detector 
feature, or to a data reduction artifact. Besides, its reality 
is confirmed by its independent detection in both the Ks- 
band images taken in the two consecutive nights as well 
as in the 7J-band image. However, the object is not de- 
tected in the J-band image and in the FORSl i?-band one 
{R ^ 25.6), where it is buried by the nebula background 
(see §2.3). Since the NACO PSF is largely oversampled, 
we have resampled all the reduced J, H and i^T^-band im- 
ages to obtain a detection with a better signal-to-noise. To 
this aim, we used the swarp program (Bertin E., Terapix 
Project) to resample the NACO images with a 3 x 3 pixels 
window while conserving the fiux. This helped in reducing 
the noise and in enhancing the faintest objects. Resampled 
images have been finally registered and co-added. Fig. [2] 
shows the co-added JHKg-hand image with the computed 
Chandra position of the CCO overlaid. The object close 
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Fig. 2. 14" X 14" section of the VLT/NACO co-added 

JHKs-h&nd image of the Vela Jr. CCO. North to the top, 
East to the left. The circle (O'.'T radius) corresponds to 
the Chandra position uncertainty after accounting for the 
accuracy of our astrometric solution (see §3.2). Star Z of 
Pavlov et al. (2001) is visible northwest of the error circle. 
The faint object {H = 21.6±0.1; K, = 21.4±0.2), labelled 
A, close to the southeast edge of the error circle is the 
possible IR counterpart of the CCO. Another faint object 
(J = 21.2 ± 0.2), labelled B, is visible ~ 1'.'9 northeast 
from the computed Chandra position. 

to the Chandra error circle is labelled A. A second point- 
like object (labelled B) is identified in the co-added image, 
although it is detected in the J-band image only, ^ 1'.'9 
northeast from the computed Chandra position. However, 
the almost 2a offset from the edge of the error circle sug- 
gests that it is likely a background object. We thus assume 
that object A is the only possible counterpart to the Vela 
Jr. CCO. 

We have computed the magnitudes of object A through 
PSF photometry. Although the accuracy of the PSF pho- 
tometry certainly depends on the AO correction, the tar- 
get faintness makes PSF photometry more suitable than 
standard aperture photometry. To this aim, we have used 
the Stetson (1992, 1994) suite of programs Daophot, fol- 
lowing the procedures outlined in Zaggia et al. (1997). 
As a reference for our photometry we have used the re- 
sampled J, H and /Cg-band images, which have a much 
better signal-to-noise. For each image we then calculated 
the model PSF with Daophot using a 6 parameters Penny 
function to fit the profile of the ^ 15 brightest non- 
saturated objects in the field. The model PSF was then 
used to measure the fluxes of object A as well as of all 
other objects in the field. To improve the precision of our 
PSF photometry, as well as to maximize the number of 
objects identified in the field, we have used the Allframe 
routine of Daophot. From the co-added JHRg-hand im- 



age we have run Daophot to create a master list of objects 
which we have properly registered on each single band im- 
age and then passed to Allframe to be used as a reference 
for PSF photometry in the single band images. The single 
band photometry catalogues have been then matched by 
Allframe to create the final multi-band photometry cat- 
alogue. Photometry measurements obtained in the two 
Ks-hand images have been averaged. We performed the 
absolute photometry calibration of our PSF photometry 
using the zero points provided by the NACO data reduc- 
tion pipeline. Since the NACO zero points are computed 
through aperture photometry, the magnitudes computed 
by Daophot have been corrected by applying the aperture 
correction. We finally applied the airmass correction us- 
ing the atmospheric extinction coefficients for the Paranal 
Observatory. The magnitudes of object A then turns out 
to be fl" = 21.6 ± 0.1 and Kg = 21.4 ± 0.2. Apart from 
object B (J = 21.2 ± 0.2), no other object is detected 
within/close to the Chandra error circle down to J ~ 22.6, 
H ~ 22.5 and K. ~ 21.8. 



4. Discussion 

4.1. The nebula 

Clumpy emission structures and filaments produced by 
shocks from the supernova explosion are detected in Ha 
in young SNRs, normally (but not only) along the external 
rims. Thus, it is possible, although very unlikely, that the 
CCO projected position is by chanche coincident with one 
of these structures. More likely, the nebula is produced by 
the neutron star ionization of the neutral Hydrogen in 
the SNR. We note that the existence of neutral Hydrogen 
in young SNRs, following atoms recombination after the 
supernova explosion, has been confirmed in a number of 
cases (see, e.g. Michael 2000). 

One of the possible interpretations of the nebula 
(Pellizzoni et al. 2002) is in terms of a bow-shock produced 
by the neutron star motion through the ISM. Interestingly, 
the i?-band image shows for the first time that the nebula 
is not spherically symmetrical but looks more like a kid- 
ney bean, a shape somehow reminescent of the arc-like 
structures observed around bow-shocks. Unfortunately, 
the spatial resolution of the VLT image is not sufficient to 
resolve more clearly the nebula morphology and to map 
its surface brightness distribution which, in the case of a 
velocity-driven bow-shock, would present a sharp enhance- 
ment towards the apex, depending on the actual value 
and direction of the neutron star velocity. We note that 
the position of the CCO lies behind the apparent apex of 
the nebula and roughly on its projected axis of simmetry, 
right where expected for a velocity-driven bow-shock. If 
object A were indeed the CCO counterpart, its position 
would be quite close to the nebula and slightly displaced 
from its apparent projected axis of symmetry. However, 
given the poor characterization of the nebula morphology 
and the lack of information on the neutron star spatial ve- 
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locity, this would not affect the vahdity of the bow shock 
interpretation. 

Since no proper motion measurement is yet available 
for the Vela Jr. CCO, the bow-shock interpretation can 
not be directly confirmed. No counterpart of the putative 
bow-shock has been detected so far in radio and in the X- 
rays. However, we note that the detection of a bow-shock 
in Ha only would not be in contradiction with its non- 
detection at other wavelengths. Indeed, the detection of a 
bow-shock in Ha is typically anti-correlated with its detec- 
tion in radio and in the X-rays, with the only noticeable 
exception of the bow-shock around the "Black Widow" 
pulsar (Stappers et al. 2003). A possible argument against 
the bow-shock interpretation comes from the position of 
the CCO with respect to the estimated geometrical center 
of the Vela Jr. SNR. According to Aschenbach (1998), this 
should be located at a( J2000)=08'* 52"^ 03", (5(J2000)= 
-46° 22', which is ~ 4 arcmin south of the CCO posi- 
tion. This would imply a proper motion direction differ- 
ent from the one expected from the apparent morphology 
of the nebula, which would be more compatible with a 
proper motion vector pointing southwest. However, both 
the uncertainty in the definition of the geometrical center 
of the SNR, see e.g. the representative case of the Vela 
SNR (Aschenbach et al. 1995), and the still poorly char- 
acterized morphology of the nebula do not make this a 
strong argument to rule out the bow-shock interpretation. 

An alternative possibility is that the nebula around 
the Vela Jr. CCO is a photo-ionization nebula, powered 
by the UV radiation from the hot neutron star. The possi- 
ble structure of a photo-ionization nebula associated to a 
neutron star moving supersonically through the interstel- 
lar medium has been studied by Van Kerkwijk & Kulkarni 
(2001). According to their model, a velocity-driven mor- 
phology is expected, with a definite symmetry with re- 
spect to the direction of the projected space velocity of 
the NS. Such a nebula should also have a smooth bright- 
ness profile, rather different from the one of a bow-shock. 
Our VLT i?-band image points to a morphology which 
could be consistent with such a picture. However, as for 
the case of the bow-shock scenario, any firm conclusion is 
prevented by the non sharp enough angular resolution of 
the ground-based images. 

Altough most of the nebula emission is in Ha (see 
§2.2), we can not a priory rule out the possibility of a syn- 
chrotron nebula, powered by the relativistic particle wind 
emitted by the neutron star. These pulsar-wind nebulae 
are often observed around young rotation-powered pul- 
sars, usually in the radio and in the X-ray bands but also in 
the optical band. Bright optical pulsar-wind nebulae have 
been indeed observed around the Crab pulsar (e.g. Hester 
et al. 1995) and PSR B0540-69 (Caraveo et al. 2001). 
However, while these objects also feature X-ray pulsar- 
wind nebulae, with structures and morphologies very simi- 
lar to the optical ones, no evidence for an X-ray nebula has 
been found for the Vela Jr. CCO. High-resolution Chandra 
X-ray observations (Pavlov et al. 2001) shows that the 
source profile is consistent with the instrument's PSF, 
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Fig. 3. (left) J vs H — Kg color magnitude diagram built 
for all objects detected in the ~ 30" x 30" field around the 
Vela Junior CCO. Field objects are marked by filled black 
squares. The possible counterpart to the CCO is marked 
by the filled red circle. The horizontal bar corresponds 
to the possible H — Kg range derived from the attached 
photometry errors. 

which rules out the presence of an extended emission com- 
ponent, unless it is extremely faint. Recent radio observa- 
tions performed with the Australian Telescope Compact 
Array (Reynoso et al. 2006) have also ruled out the pres- 
ence of a pulsar-wind nebula in radio. Thus, the lack of an 
X-ray/radio counterpart, together with the apparent lack 
of continuum optical emission, makes the pulsar-wind neb- 
ula interpretation for the Vela Jr. nebulosity very unlikely. 

4.2. The point source 

At the moment we can not confidently rule out the possi- 
bility that object A is just a field object unrelated to the 
CCO, whose IR counterpart would thus remain unidenti- 
fied. Indeed, given the density of stars in the NACO field 
of view we estimate a chance coincidence probability of 
~ 15%, i.e. certainly not negligible. Given our flux upper 
limits {R ~ 25.6, J - 22.6, H - 22.5 and Ks ~ 21.8), 
this would almost certainly rule out an hypothetical stel- 
lar companion (see discussion below) and would instead 
support the conclusion that the CCO is an isolated neu- 
tron star, possibly surrounded by a circumstellar disk. 
However, if object A is indeed the CCO counterpart, its 
nature is controversial. As the matter of fact, it is not 
a priori obvious whether the detected IR emission can 
be ascribed to the neutron star itself, to a fallback disk 
around the neutron star, or to a putative companion star. 
Suggestive as it could be, the object does not show evi- 
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dence for variability, with the Kg magnitudes measured in 
the two nights consistent within the associated errors. 

The first possibility is that object A is itself the IR 
counterpart of the neutron star. So far, IR emission has 
been detected only for five rotation-powered pulsars (see 
Mignani et al. 2007b for a summary) and it is ascribed to 
synchrotron radiation produced by relativistic particles in 
the neutron star's magnetosphere. The observed IR lumi- 
nosity spans about six decades and it is strongly depen- 
dent on the neutron star's age. If we assume for the Vela 
Jr. CCO an age of a few thousands years, as estimated 
for the host SNR, it might then have an IR luminosity 
comparable to the one of the Vela pulsar, i.e. ~ 10^^ ergs 
s~^. At a distance of 1 kpc it would then have a magni- 
tude Kg ~ 24.1, i.e. much fainter than the one of object 
A. IR emission from an isolated neutron would thus be in- 
compatible with our detections unless the CCO distance 
is as low as ~ 300 pc, or its IR luminosity is a factor ~ 10 
higher than expected. However, given both the uncertain- 
ties on the CCO distance and the lack of information on 
its energetics, we can not a priori exclude that object A is 
itself the counterpart of the neutron star. 

An alternative possibility is that the IR emission from 
object A is produced from a fallback disk around the neu- 
tron star. Fallback disks have also been invoked to ex- 
plain the IR spectral flattening of the AXPs (e.g., Israel 
et al. 2003; Israel et al. 2004) and have been searched 
around other objects possibly related to the magnetars, 
like the high-magnetic field radio pulsars (Mignani et al. 
2007b) and the XDINSs (Lo Curto et al. 2007) but, so far, 
with unconclusive results. Evidence for a fallback disk has 
been possibly found from Spitzer observations of the AXP 
4U 0142-H61 (Wang et al. 2006). We note that the Fir/ Fx 
ratio of the Vela Jr. CCO is ~ 5 x 10^^, which is very simi- 
lar to the value observed for the case of AXP 4U 0142-1-61. 
Accurate multiband IR photometry is required in order to 
test such an hypothesis. 

The last possibility is that object A is the stellar mass 
companion of the CCO. We have tried to determine a pos- 
sible stellar classification of object A by using the available 
color information. Unfortunately, its non-detection in the 
J band leaves us with only a measure oi H — Kg and with 
upper limits on J — H and J ~ Kg. However, the fact that 
object A has not been detected in the J band but only 
in the H and Kg bands, despite of similarly deep expo- 
sures, seem to imply that its spectrum is quite red. To 
better evaluate its colors, we have compared its location 
in a J vs. ff — Kg color magnitude diagram (CMD) with 
respect to those of a number of stars detected in the field 
(Fig. [3K). Photometry of the field stars has been com- 
puted as described in §3.3. Unfortunately, only ~ 50 stars 
have been identified in the narrow and sparsely populated 
NACO field of view, which certainly hampers our char- 
acterization of the field star sequence. At a first glance, 
with an H — Kg ^0.2, object A indeed seems to devi- 
ate from the sequence. However, its large H — Kg error 
makes this deviation not significant. Thus, object A is not 
redder, possibly only slightly bluer, than the majority of 



the field stars. This can be visually appreciated from the 
composite JHKg-hand image of the field (Fig. [Sb). We 
have investigated the effects of the interstellar extinction 
on the CMD. The interstellar reddening towards the Vela 
Jr. CCO can be estimated from the Hydrogen column den- 
sity A^H ~ 3.5 10^^ cm— 2 derived from the X-ray observa- 
tions (Kargaltsev et al. 2002). According to the relation of 
Predhel & Schmitt (1995) this gives an E(B - F) ~ 0.67 
which corresponds to Aj ~ 0.5, Ah ~ 0.3 and Ak ~ 0.2 
assuming the extinction coefficients of Fitzpatrick (1999). 
However, accounting for the reddening does not affect sig- 
nificantly the location of object A in the CMD with re- 
spect to the sequence of the field stars. Thus, we conclude 
that object A has no peculiar colors and it does not stand 
out from the average stellar population in the field. 

The observed colors of object A, if - ivTs = 0.2 ± 0.2, 
J — H > 1.0, and R — Kg > 5.2 can be compatible with 
those of a mid M type star, after accounting for the in- 
terstellar reddeding (e.g. Legget 1992). For instance, an 
M5 star at a distance of 2 kpc would be just below our 
detection limit in J but would have H and i^s-band mag- 
nitudes consistent with the ones measured for object A. 
We note that a distance of 2 kpc is larger than the one 
of 1 kpc usually assumed for the SNR (Slane et al. 2001), 
but is not incompatible with the observations. Thus, our 
data are consistent with the Vela Jr. CCO being a neu- 
tron star in a binary system with an M type companion 
star. However, since we have only a lower limit on J — H , 
we can not exclude that object A is actually a star of an 
even later spectral type, an example of the recently dis- 
covered L dwarfs (Kirkpatrick et al. 1999). Late type star 
counterparts have been also proposed in the past for the 
PKS 1209-52 (Pavlov et al. 2004) and RCW103 (Sanwal 
et al. 2002) CCOs but both identifications have been later 
discarded from a re-analysis of the Chandra astrometry 
(Mignani et al. 2007a; Wang et al. 2007; De Luca et al. 
private communication). If our identification is confirmed, 
the association of a young neutron star in a binary system 
with such a low mass companion will pose more than one 
theoretical challenge for current neutron star formation 
and supernova explosion models. 

5. Conclusions 

We have reported on optical/IR VLT observations of the 
Vela Jr. CCO. The VLT i?-band data have revealed a com- 
pact nebula at the Chandra position of the CCO, consis- 
tent both in morphology and in brightness with the Ha 
nebula identified by Pellizzoni et al. (2002). We can not 
presently determine whether the nebula is associated to a 
bow-shock produced by the supersonic motion of the neu- 
tron star in the ISM, or it is a photo-ionization nebula 
produced by the UV thermal radiation from the neutron 
star. Higher resolution imaging is required to resolve the 
nebula structure and to determine whether it is more com- 
patible with the bow-shock or the photo-ionization sce- 
nario, for which we expect quite distinct morphologies, 
i.e. arc-like or spherical, respectively. Alternatively, high- 
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resolution spectroscopy would make it possible to mea- 
sure the width of the hue and to derive its velocity 
broadening which is typically higher for a velocity driven 
bow-shock nebula. 

Proper motion measurement of the CCO also represent 
a powerful diagnostic to investigate the bow-shock sce- 
nario, where we expect the velocity vector to lie along the 
simmetry axis of the nebula and to point towards its apex. 
Since the CCO is radio quiet, one possibility to obtain a 
proper motion measurement would be through high reso- 
lution X-ray observations with Chandra. However, proper 
motion measurement performed with Chandra on some 
nearby INSs (Motch et al. 2007) only allowed to obtain 
marginally significant results. For the Vela Jr. CCO which 
is up to a factor five more distant, such measurements are 
even more difficult. More realistically, a proper motion 
measurement could be obtained from high-resolution IR 
imaging of the candidate CCO counterpart identified in 
our VLT H and Kg-hand observations. Its proper motion 
measurement will allow to test the bow-shock interpreta- 
tion for the nebula on one hand, and to confirm its pro- 
posed identification, which only relies on the coincidence 
with the Chandra position, on the other. 

The nature of the candidate CCO counterpart is, at 
present, unclear. It could be the neutron star itself, or a 
fallback disk left over from the supernova explosion, or a 
very late type star in a binary system with the neutron 
star. In all cases, confirming the identification will have 
several important implications for our understanding of 
the CCOs as class, as well as for neutron star formation 
and supernova explosion models. Instead, discarding the 
proposed identification will strengthen the more standard 
framework isolated neutron star scenario for the CCOs. 
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